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Abstract 

A new class of the constrained Monte Carlo (CMC) algorithms for the QCD 
evolution equation was recently discovered. The constraint is imposed on the type 
and the total longitudinal energy of the parton exiting QCD evolution and entering 
a hard process. The efficiency of the new CMCs is found to be reasonable. 
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1 Introduction 



This brief report summarizes the recent developments in the area of the Monte Carlo 
(MC) techniques for the perturbative QCD calculations. Most of it was done at the time 
of the present HERA-LUC workshop, partial results being presented at several of its 
meetings. At present, two papers, £Q and j2], demonstrating the principal results are 
already available. Generally, these MC techniques concern the QCD evolution of the 
parton distribution functions (PDFs) D k (x,Q), where k denotes the type of the parton 
(quark, gluon), x the fraction of longitudinal momentum of the initial hadron carried 
by the parton, and the size of the available real/virtual emission phase space is Q. The 
evolution equation describes the response of the PDF to an increase of Q; Dk{x, Q) is an 
inclusive distribution and can be measured almost directly in hadron-lepton scattering. 
On the other hand, it was always known that there exists in QCD an exclusive picture 
of the PDF, the so-called parton-shower process, in which Dk(x,Q) is the distribution 
of the parton exiting the emission chain and entering the hard process (lepton-quark for 
example). The kernel functions Pkj(Q, z), that govern the differential evolution equations 
of PDFs are closely related to distributions governing a single emission process (i — 1) — ► i 
in the parton shower: Pk^AQi^H 

In other words, the evolution (Q- dependence) of PDFs and the parton shower represent 
two faces of the same QCD reality. The first one (inclusive) is well suited for basic precision 
tests of QCD at hadron-lepton colliders, while the second one (exclusive) provides realistic 
exclusive Monte Carlo modelling, vitally needed for experiments at high-energy particle 
colliders. 

At this point, it is worth stressing that, so far, we were referring to DGLAP-type 
PDFs [3] and their evolution, and to constructing a parton-shower MC starting from 
them, as was done two decades ago and is still done today. This involves a certain 
amount of "backward engineering" and educated guesses, because the classical inclusive 
PDFs integrate over the pt of the exiting parton. The so-called unintegrated PDFs 
(UPDFs) Dk(x,pT,Q) would be more suitable for the purpose, leading to higher- quality 
QCD calculations. UPDFs are, however, more complicated to handle, both numerically 
and theoretically. (It is still a challenge to construct a parton-shower MC based consis- 
tently on the theoretically well defined UPDFs.) 

Another interesting "entanglement" of the evolution of PDFs on one side and of the 
parton shower (PS) MC on the other side is also present in the modelling of the showering 
of the incoming hadron - mostly for technical reasons and convenience. The Markovian 
nature of the QCD evolution can be exploited directly in the PS MC, where partons 
split / decay as long as there is enough energy to dissipate (final state) or the upper bound- 
ary Q of the phase space is hit (initial state). The multiparton distribution in such a MC 
is a product of the evolution kernels. However, such a direct Markovian MC simulation 
of a shower is hopelessly inefficient in the initial state, because the hard process accepts 
only certain types and momenta of the incoming partons - most of the shower histories 
are rejected (zero MC weight) by the hard process. This is especially true for forming 
narrow resonances such as electroweak bosons or Higgs boson at the LHC. The well known 
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"workaround" is the backward evolution MC algorithm of Sjdstrand, used currently in all 
PS MCs, such as HERWIG [4 and PYTHIA 0. Contrary to the forward Markovian MC, 
where the physics inputs are PDFs at low Qq ~1 GeV and the evolution kernels, in the 
backward evolution MC one has to know PDFs in the entire range (Qo, Q) from a separate 
non-MC numerical program solving the evolution equation and setting up look-up tables 
(or numerical parametrization) for them 1 . 

The following question has been pending in the parton-shower MC methodology for a 
long time: Could one invent an efficient "monolithic" MC algorithm for the parton shower 
from the incoming hadron, in which no external PDFs are needed and the only input are 
PDFs at Qo an d the evolution kernel (the QCD evolution being a built-in feature of the 
parton shower MC)? Another question rises immediately: Why bother? Especially since 
this is a tough technical problem. This cannot still be fully answered before the above 
technique is applied in the full-scale (four-momentum level) PS MC. Generally, we hope 
that this technique will open new avenues in the development of the PS MC at the next- 
to-leading-logarithmic (NLL) level. In particular, it may help in constructing PS MCs 
closely related to unintegrated structure functions and, secondly, it may provide a better 
integration of the NLL parton shower (yet to be implemented!) with the NLL calculation 
for the hard process. 

The first solution of the above problem of finding an efficient "constrained MC" (CMC) 
algorithm for the QCD evolution was presented in refs. PUE1- This solution belongs to 
what we call a CMC class II, and it relies on the observation that all initial PDFs at Qq 
can be approximated by const • Xq _1 ; this is to be corrected by the MC weight at a later 
stage. This allows elimination of the constraint x = Yii z ii at the expense of xq, keeping 
the factorized form of the products of the kernels. Simplifying phase-space boundaries 
in the space of Zi is the next ingredient of the algorithm. Finally, in order to reach a 
reasonable MC efficiency for the pure bremsstrahlung case out of the gluon emission line, 
one has to generate a 1/z singularity in the G —>■ G kernel in a separate branch of the 
MC. The overall efficiency of the MC is satisfactory, as is demonstrated in ref. pQ for the 
case of the pure bremsstrahlung out of the gluon and quark colour charge. Generalization 
to the quark-gluon transition is outlined, but not yet implemented. The main drawback 
of this method is its algebraic complexity. Further improvement of its relatively low MC 
efficiency is possible (even though it could lead to even more algebraic complexity). 

The second, more efficient, CMC algorithm was presented in ref. [2] (as well as during 
the October 2004 meeting of the workshop). It belongs to what we call a CMC class I. 
The main idea is to project/map points from the hyperspace defined by the energy con- 
straint x = rj. Zi, into a simpler hyperspace, defined by the hardest emission, x = minZj. 

backward evolution is basically a change in the order of the generation of the variables: Consider 
generating p(x,y), where one generates first x according to p{x) = J dy p{x,y), and next y according 
to p(x,y), by means of analytical mappings of x and y into uniform random numbers. However, such 
analytical mappings may not exist, if we insist on generating first x and next y\ Nevertheless, we may still 
proceed with the same method by "brute force" , if we prctabulatc and invert numerically the functions 
R{x) = J x J dx'dy' p(x', y') and R x (y) = J V dy' p(x, y'). This is what is done in a more dimensional case 
of the backward-evolution MC; it also explains why pretabulated PDFs are needed in these methods. 
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This mapping is accompanied by the appropriate MC weight, which compensates exactly 
for the deformation of the distributions involved, and the bookkeeping of the hyperspace 
boundaries is rigorous. The above describes a CMC for the pure bremsstrahlung seg- 
ment of the gluon emission out of a quark or gluon chain. Many such segments are 
interconnected by the quark-gluon transitions. The algebraic hierarchic reorganization 
of the emission chain into a super-level of the quark-gluon transitions and sub-level of 
the pure bremsstrahlung is an important ingredient in all CMC algorithms and will be 
published separately [ZJ . The basic observation made in ref . [Hj is that the average number 
of super-level transitions is low, ~ 1; hence for precision of a 10~ 4 it is sufficient to limit 
it to three or four transitions. The integration/simulation of the super-level variables is 
done efficiently using the general-purpose MC tool FOAM El] . The above proof of 
the correctness of the CMC class I algorithm concept was given in ref. [2] for the full 
DGLAP-type QCD evolution with the LL kernels (including quark-gluon transitions). 

Although our main aim is to construct the non-Markovian CMC class of algorithms, 
we have developed in parallel the family of Markovian MC (MMC) algorithms/programs, 
which provide numerical solutions of the QCD evolution equations with high precision, 
~ 10~ 3 . We use them at each step of the CMC development as numerical benchmarks 
for the precision tests of the algorithms and their software implementations. The first 
example of MMC for DGLAP at LL was defined/examined in ref. jH] and tested using the 
non-MC program QCDnuml6 11 2 . In some cases our MMC programs stand ahead of 
their CMC brothers; for instance, they already include NLL DGLAP kernels. A systematic 
description of the MMC family of our MC toolbox is still under preparation [13J. 

The last development at the time of the workshop was an extension of the CMC 
type-I algorithm from DGLAP to CCFM one-loop evolution (also referred to as 
HERWIG evolution in which the strong coupling constant gains z-dependence, 

ct s {Q) —> o( s (Q(l — z)), as advocated in ref. [TH], confirmed by NLL calculations [T7] . 
The above ansatz also compels introduction of a Q-dependent IR cutoff, e = Q e /Q: an- 
other departure from DGLAP. This version of the CMC is still unpublished. Its version 
for the pure bremsstrahlung was presented at the March 2005 meeting of the workshop; in 
particular a perfect numerical agreement with the couterpartner MMC was demonstrated. 
Recently both CMC and MMC for the one-loop CCFM were extended to quark-gluon 
transitions, and again perfect agreement was found. 

For the detailed description of the new CMC algorithm, we refer the reader to the 
corresponding papers pQ and [2j and workshop presentations 3 . Here, let us only show 
one essential step in the development of the CMC for the one-loop CCFM model - the 
mapping of the Sudakov variables for the pure bremsstrahlung: 

dt / dz a(Q(l - z)) zP^ G (z, t) 

= i r 1 f * = I r *« f *(«)■ 

ft Jo A, t + ln( l -z) l-z I3j n ' J„ w 

2 It was also compared with the non-MC program APCheb |12) . 
3 To be found at http : // j adach . home . cern . ch/ j adach/ 
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Figure 1: CMC of the one-loop CCFM versus the corresponding MMC for quarks; number of 
quark-gluon transitions J = 0,1,2,3,4, and the total. The ratios in the lower plot are for 
n = 0, 1 and the total (blue). 

The short-hand notation i = i(t) =t — t^ and v = ln(l — z) supplements that of ref. [2j 
in use, and the mapping reads 

y{z) = p[vi;ti,t ) = p(yi +t x ) - 9 Vl>te - to p{vi +t ), s(t) 



p'(v-i u i o y (2) 
p'(v ; ti, i ) = 9 v<tE ^ to p'(v + t{) + 9 v>te -t [p'{v + h) - p'{v + i )), 

where p(t) = t(lnt — lnt e ) + 1 £ — t. Once the above mapping is set, the same algo- 
rithm, with the parallel shift i/i — > y, t + Y, can be used in this case. The super-level of 
quark-gluon transitions is again implemented using FOAM 4 . A numerical comparison of 
the corresponding CMC and MMC programs is shown in fig. The MC efficiency is 
comparable with that of the DGLAP case. 

Summary: We have constructed and tested new, efficient, constrained MC algorithms 
for the initial-state parton-emission process in QCD. 
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